Journal  of  Power  Sources  195  (2010)  1383-1392 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


in 

SllbaidtS 


A  high  efficient  assembly  technique  for  large  proton  exchange  membrane  fuel 
cell  stacks:  Part  II.  Applications 

P.  Lin,  P.  Zhou,  C.W.  Wu* 

State  Key  Lab  of  Structural  Analysis  for  Industrial  Equipment,  Faculty  of  Vehicle  Engineering  and  Mechanics,  Dalian  University  of  Technology,  Dalian,  Liaoning  1 1 6024,  China 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  15  September  2009 
Accepted  15  September  2009 
Available  online  28  September  2009 


Keywords: 

Fuel  cell 

Proton  exchange  membrane 
Clamping  load 
Assembly  technique 


A  high  efficient  assembly  technique  for  large  proton  exchange  membrane  fuel  cell  (PEMFC)  stacks  is 
proposed  for  obtaining  the  optimal  clamping  load.  Using  the  equivalent  stiffness  model  proposed  in  Part 
I  of  this  study,  we  show  how  to  design  the  structure  components  for  a  large  PEMFC  stack.  First,  we  give 
a  design  demonstration  based  on  the  structural  strength  of  the  stack.  We  then  discuss  how  to  obtain  the 
optimal  clamping  load  for  a  given  PEMFC  stack  according  to  the  requirements  of  the  interface  contact 
resistance  and  permeability  of  the  gas  diffusion  layer.  Finally,  we  discuss  the  effects  of  the  equivalent 
stiffness  of  the  spring  washer  on  the  structure  thermal  stress. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  have  many 
advantages  such  as  their  clean  and  quiet  operation,  high  efficiency, 
low  operation  temperature,  and  quick  start-up,  and  they  are  there¬ 
fore  a  possible  alternative  to  combustion  engines.  Studies  on  a 
single  PEMFC  during  the  past  decade  have  given  a  very  good  under¬ 
standing  of  its  basic  principles.  However,  in  the  last  few  years, 
the  performance  of  a  large  PEMFC  stack  has  attracted  more  and 
more  attention.  To  design  a  fuel  cell  stack  with  high  reliability  and 
long  lifetime,  there  are  a  number  of  new  technologies  that  need 
to  be  developed.  Among  them,  the  assembly  is  known  to  be  a  key 
technique,  with  previous  studies  having  shown  that  the  clamping 
load  (clamping  pressure)  of  a  single  PEMFC  plays  an  important  role 
in  optimizing  its  performance  [1-7].  In  Part  I  of  this  study  [8],  a 
high  efficient  assembly  model  was  proposed,  where  a  large  PEMFC 
stack  system  was  considered  as  a  mechanical  equivalent  stiffness 
model  consisting  of  numerous  elastic  elements  (springs)  connected 
in  either  series  or  parallel.  Comparison  with  a  three-dimensional 
(3D)  finite  element  analysis  (FEA)  showed  that  the  equivalent  stiff¬ 
ness  model  had  very  good  prediction  accuracy  for  calculating  the 
component  stiffness  and  clamping  load.  In  the  present  paper,  prac¬ 
tical  examples  are  introduced  to  show  how  to  design  the  optimal 
clamping  load  for  a  large  PEMFC  stack  using  the  equivalent  stiff¬ 
ness  model.  However,  the  basic  idea  can  be  used  in  designing  other 
types  of  fuel  cells  as  well. 
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The  clamping  load  of  a  large  PEMFC  stack  affects  the  perfor¬ 
mance  and  lifetime  in  several  ways.  Too  large  a  clamping  load 
may  cause  some  components  in  the  stack  to  produce  a  stress  high 
enough  to  give  rise  to  plastic  deformation  or  even  crack  [9],  while 
too  small  a  clamping  load  may  result  in  a  high  contact  electrical 
resistance  at  the  interface  between  the  gas  diffusion  layer  (GDL) 
and  the  bipolar  plate  (BPP)  [1-3]  or  the  leakage  of  water  or  fuel 
at  the  sealed  interfaces  [10].  Even  if  a  reasonable  initial  clamping 
load  is  applied  at  room  temperature,  either  the  contact  electrical 
resistance  or  the  sealed  interface  pressure  or  the  structural  stress 
may  fail  to  meet  the  design  criteria  when  considering  the  thermal 
effects  under  working  conditions.  Therefore,  the  thermal  effects 
should  be  taken  into  consideration  when  calculating  the  pressure 
distribution  inside  PEMFC  stacks.  However,  the  computing  time  for 
such  a  large  stack  using  a  numerical  method  (usually  FEA)  is  usu¬ 
ally  unacceptable  owing  to  the  complexity  of  the  stack  structure.  In 
the  past,  a  few  attempts  have  been  made  to  simulate  the  compres¬ 
sion  pressure  within  a  single  PEMFC  [4,1 1-13].  Even  for  a  single  cell 
computational  model,  it  usually  takes  a  few  hours  for  a  step  com¬ 
putation.  In  the  following,  we  give  a  detailed  description  of  how 
to  obtain  solutions  for  the  problems  mentioned  above  using  the 
spring  assembly  model  proposed  in  the  previous  study  [8]. 

2.  Pretension  and  tightening  torque  of  clamping  bolts 

In  the  traditional  PEMFC  stacks,  end  plates  (EPs)  together  with 
an  array  of  bolts  are  usually  used  to  fasten  the  inner  components, 
as  shown  in  Fig.  1.  This  connecting  method  is  called  the  “tie-rod 
and  bolt”  method  [14].  The  clamping  load  is  controlled  through 
adjusting  the  tightening  torque  applied  to  the  nuts. 
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Nomenclature 

A  cross-section  area  (m2) 

b  width  of  the  sealant  (m) 

C  design  number  of  bolts  in  the  stack 

£h  outside  diameter  of  the  nut  (m) 

D0  inside  diameter  of  the  nut  (m) 
d  nominal  diameter  of  the  bolt  (m) 

d\  minor  diameter  of  the  bolt  (m) 

deq  equivalent  frictional  diameter  (m) 

F  load  (N) 

k  equivalent  stiffness  (N  nrr1 ) 

M  design  number  of  cells  in  the  stack 

nsafe  safety  factor 

P  pressure  (Pa) 

T  tightening  torque  (N  m) 

Z thermai  thermal  deformation  (m) 

Greek  symbols 

P  half  of  the  thread  angle  in  the  vertical  cross-section 

(°) 

y  number  of  ribs  in  a  single  flow  field 

8  thread  pitch  (m) 

£  coefficient  of  material  transience  from  elastic  to 

plastic 

0  frictional  angle  (°) 

k  in-plane  gas  permeability  of  the  GDL  (m2) 

/x  friction  coefficient  between  a  screw  pair 

p  contact  resistivity  of  the  GDL  ( £2  m2 ) 

a  stress  (Pa) 

cp  slope  angle  of  the  thread  (°) 


Fig.  2.  Schematic  diagram  of  the  forces  acting  on  the  rectangular  thread:  (a)  the 
rectangular  thread  and  (b)  the  forces  acting  on  the  rectangular  thread. 


When  tightening  a  bolt  using  a  wrench,  the  total  torque  applied 
to  the  nut  equals  the  sum  of  the  frictional  torque  between  the  screw 
pair  (Ti )  and  that  between  the  nut  and  EP  surface  (T2): 

T  =  T\  +T2.  (1) 


where  <p  is  the  slope  angle  of  the  thread,  F \  is  the  necessary  force 
for  spiral  ascending,  F2  is  the  pretightening  force,  /x  is  the  friction 
coefficient  between  the  screw  pair,  and  8  is  the  thread  pitch.  Letting 
6  be  the  frictional  angle  corresponding  to  /x  (/x  =  tan0),  Eq.  (2)  can 
be  rewritten  as 


For  the  rectangular  thread  structure  shown  in  Fig.  2,  we  have 


F!  =  F2  tan(<p  +  0). 


(3) 


Fi  cos  (p  -  F2  sin  (p  -  /x(Fi  sin  (p  +  F2  cos  (p)  =  0, 


(2) 


Thus,  the  frictional  torque  applied  to  the  screw  pair  is 

Fid  F2d  tan(<p  +  0) 

T,  =  —  = - 2 - ’ 


where  d  is  the  nominal  diameter  of  the  bolt. 

On  the  other  hand,  the  frictional  torque  between  the  nut  and  EP 
surface  is 


T2 


MnF2  rD' 

n{D\-Dl)JDQ 


7tD2dD  = 


3(D2-D2)  ’ 


(5) 


where  D !  is  the  outer  diameter  of  the  nut,  D0  is  the  inner  diameter  of 
the  nut,  and  /xn  is  the  coefficient  of  the  friction  between  the  contact 
surfaces  (see  Fig.  3).  Therefore,  the  relation  between  the  pretension 
force  and  tightening  torque  is 


T  =  h  +  T2  =  F2 


d  tan(y  +  9)  /MDp  -  Dp3) 
2  3(D,2-D02) 


(6) 


In  general,  the  magnitudes  of  cp  and  0  are  extremely  small.  In 
such  a  case,  we  have 


tan(<p  +  0)  «  tan  <p  +  tan  6.  (7) 

In  addition,  the  thread  pitch  can  be  expressed  as 
8  =  nd  tan  (p. 


(3) 
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Fig.  3.  Schematic  diagram  of  the  nut  geometry. 


Thus,  Eq.  (6)  can  be  rewritten  in  an  approximate  form  as 

T^F  [  8  /id  nn(Dp-D03) 

~  2  L2*  2  +  3 (D!2-D02)  _  ‘ 


O) 


For  the  triangle  thread  structure,  Eq.  (9)  should  be  modified  as 


T^F2 


8  fid  fin(D]-D30) 
27T  +  2  cos^  +  3 (D\  -Dg) 


(10) 


where  is  half  of  the  thread  angle.  In  the  following,  only  triangle 
thread  structures  are  adopted.  In  addition,  a  parameter  called  the 
“equivalent  frictional  diameter”  (deq)  is  introduced  here  for  equa¬ 
tion  simplification: 


,  _  <$  fid  fin(D3-D3) 

eq~2^  +  2~coi^0  +  3(D2-D2)  ' 


(11) 


Therefore,  Eq.  (10)  can  be  rewritten  for  PEMFC  stacks: 


T 


Fclampingdeq 

c 


(12) 


where  Fdamping  is  the  summation  of  the  clamping  loads  applied  to 
different  bolts  and  C  is  the  total  number  of  tightening  bolts  in  the 
stack. 


3.  Stack  design  based  on  structural  strength 

As  the  voltage  and  power  generated  by  a  single  PEMFC  are  rather 
limited,  fuel  cell  products  in  engineering  always  comprise  numer¬ 


Fig.  5.  Schematic  diagram  of  the  bolt  connection  using  washers  and  nut. 


ous  single  cells  in  series  that  form  a  large  stack  with  a  cyclic  cell 
characteristic  as  shown  in  Fig.  4.  Although  the  working  principle  is 
not  different  for  a  single  cell  and  a  stack,  influences  of  the  clamping 
load  and  thermal  effects  on  the  system  performance  are  more  com¬ 
plex  for  the  stack  than  for  the  single  cell  [13].  In  the  following,  we 
describe  in  detail  the  design  procedure  for  a  demonstration  PEMFC 
stack  structure. 

A  fuel  cell  stack  generally  comprises  two  types  of  structures:  (a) 
structures  used  to  generate  power,  usually  residing  in  the  inter¬ 
nal  region  of  the  stack,  and  (b)  structures  for  the  sealing  function, 
which  are  usually  located  in  the  external  region  of  the  stack.  We 
refer  to  them  as  the  internal-stack  and  external-stack  respectively 
in  the  equivalent  stiffness  model  [8].  For  the  internal-stack,  adjust¬ 
ing  the  clamping  load  usually  influences  power-related  parameters 
such  as  the  permeability  and  interfacial  resistivity  of  the  GDL.  For 
the  external-stack,  however,  adjusting  the  clamping  load  mainly 
influences  the  sealing  performance. 

The  following  design  description  is  mainly  based  on  the  struc¬ 
tural  strength.  Unless  noted  otherwise,  all  the  symbols  used  in 
this  paper  have  the  same  meanings  as  those  used  in  our  previous 
report  [8]. 

3.1.  Relation  between  the  tightening  torque  and  clamping  load 

Washers  are  a  type  of  fastener  usually  used  with  bolts  and  nuts. 
They  are  generally  divided  into  flat  washers  and  spring  washers, 
as  shown  in  Fig.  5.  If  washers,  especially  spring  washers,  need  to 
be  used  in  a  stack  assembly,  the  original  equivalent  stiffness  model 
established  in  our  previous  report  [8]  may  need  a  simple  amend- 


End  Plate 

Bipolar  Plate 

y  f  Sealants 

^  J  Membrane  Electrode  Assembly 


Clamping  Force 


Fig.  4.  Schematic  diagram  of  the  PEMFC  stack  hardware. 
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C  x  kwasher 
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End  Plate 


ts  /\ 


9  k internal-stack  k external-stack 


Tensile 

direction 

C  ^  kbolt 


End  Plate 


/////////// 


Fig.  6.  Schematic  diagram  of  the  modified  equivalent  stiffness  model  when  con¬ 
sidering  the  washer  deformation.  (The  spring  dimension  does  not  represent  the 
magnitude  of  stiffness.) 


ment,  as  shown  in  Fig.  6.  Moreover,  the  elastic  deformation  of  nuts 
is  not  considered  in  the  equivalent  stiffness  model  owing  to  their 
high  stiffness  and  limited  thickness. 

According  to  the  rigid  EP  assumption,  the  internal  and  external- 
stacks  have  the  same  deformation  after  applying  the  clamping  load; 
that  is, 

F internal-stack  _  F external-stack  _  _ F whole-stack _ 

^internal-stack  ^external-stack  ^internal-stack  +  ^external-stack 

(13) 

From  the  conditions  of  static  load  equilibrium,  one  obtains 

F whole-stack  =  F damping  •  (14) 

The  left  side  of  Eq.  (14)  represents  the  total  compression  load 
applied  to  the  stack,  while  the  right  side  is  the  total  tensile  load 
applied  to  all  bolts.  Those  two  loads  are  opposite  in  direction,  but 
the  same  in  magnitude.  Substituting  Eqs.  (14)  into  (13),  one  has 


clamping 


(^internal- stack  ~t~  ^external-stack )F internal-stack 
^internal-stack 

( ^internal-stack  +  kextemal-stack)Fexternal-stack 
^external-stack 


(15) 


Combining  Eq.  (12)  with  (15)  gives  relations  between  the  tight¬ 
ening  torque  and  clamping  load: 


T  ^  (^internal-stack  +  ^external-stack)F internal-stack  j 

1  ~  - - Oeq 

internal-stack 

and 


T  ^  (^internal-stack  +  ^external-stack  )F external-stack  j 
1  ~  7T,  ueq. 

external -stack 


(16) 


(17) 


3.2.  Clamping  load  design  based  on  the  strength  of  the 
external-stack 


For  practical  PEMFCs,  sealants  (e.g.,  sealing  washers  or  sealing 
rings  [15])  are  usually  placed  between  BPPs  to  prevent  the  leak¬ 
ing  of  gas  and  coolant,  the  typical  materials  of  which  are  Viton, 
ethylene-propylene-diene  monomer  and  silicone  [16].  One  of  the 
most  important  factors  affecting  the  sealing  performance  is  the 
mean  contact  pressure  on  the  sealing  surfaces.  Too  low  a  con¬ 
tact  pressure  may  cause  leakage,  while  too  high  a  contact  pressure 
may  damage  the  sealants  or  other  components.  To  ensure  sealing 
performance,  the  mean  contact  pressure  ( Pseaiant )  should  meet  the 


following  condition. 

pmin  ^  n  pmax 

1  sealant  —  1  sealant  —  1  sealant ’ 


(18) 


where  P^ant  is  the  necessary  mean  contact  pressure  to  prevent 
leakage  and  P™%jant  is  the  maximum  mean  contact  pressure  for 
which  the  sealant  material  does  not  fail.  From  experimental  studies 
[17],  we  have  the  following  empirical  equations. 


pmin 

1  sealant  ~ 


5.5  X  106  +  3.2Pgas 
VlOOOb 


(19) 


pmax  _  foyieM 
1  sealant  *>u  sealant  ’ 


(20) 


where  Pgas  is  the  gas  pressure  of  the  fuel  cell,  b  is  the  width  of  the 
sealant,  £  is  the  transient  coefficient  of  the  sealant  material  from 
elastic  to  plastic  deformation  (roughly  having  a  value  of  3  for  most 
polymeric  materials),  and  is  the  yield  stress  of  the  sealant 

material.  The  range  of  the  clamping  load  applied  to  the  external- 
stack  can  then  be  determined  using 


(5.5  x  10  +  3 .2Pgas)Aseaiant 

V1000 b 


yyield 

J  sealant1 


(21) 


Therefore,  the  allowed  range  of  the  tightening  torque  applied  to  a 
single  bolt  on  the  basis  of  the  strength  requirement  of  the  external- 
stack  can  be  written  as 


T  > 


(5.5  X  10  +  3.2Pgas)(kjnferna/_stack  T  kexternal-StaCk)Asealantdeq 


V1000 bCkextemai 

-stack 


T  < 


(^internal-stack  +  ^externa/-stack)^°rJea/anf71sea/ant^eq 
Ckextemal-stack 


3.3.  Clamping  load  design  based  on  the  strength  of  the 
internal-stack 


(22) 

(23) 


The  membrane  electrode  assembly  (MEA)  plays  an  important 
role  in  power  output.  It  is  located  in  the  internal-stack  and  can  be 
considered  a  sandwich  structure  of  “GDL  + proton  exchange  mem¬ 
brane  (PEM)  +  GDL”  [8].  Since  the  GDL  and  PEM  can  be  both  highly 
deformed  owing  to  either  the  high  porosity  or  the  low  elastic  mod¬ 
ulus,  the  MEA  deformation  makes  a  primary  contribution  to  the 
stack  deformation  [18].  Most  of  the  MEA  deformation  is  caused  by 
the  clamping  load,  although  other  factors,  such  as  the  pressures  and 
humidities  of  the  fuel  and  oxidant,  the  material  properties,  and  the 
thickness  of  the  sealant  layer,  may  also  affect  the  MEA  deformation. 
It  has  been  reported  that  failures  of  MEAs  are  mainly  caused  by  the 
clamping  load  and  thermal  stress  [19,20]. 

An  electric  circuit  forms  as  soon  as  the  BPP  makes  contact  with 
the  GDL  during  assembly.  If  the  clamping  load  applied  to  the  contact 
surface  is  not  large  enough,  there  is  an  unreasonably  high  contact 
resistance  that  produces  an  unexpected  voltage  drop  in  the  circuit. 
Therefore,  the  minimum  clamping  load  applied  to  the  internal- 
stack  should  be  higher  than  the  load  caused  by  the  fuel  gas  pressure. 
Elowever,  it  can  share  the  same  lower  bound  with  the  external- 
stack  as  given  by  Eq.  (22).  On  the  other  hand,  the  clamping  load 
applied  to  the  internal-stack  should  not  be  so  large  that  the  stress 
in  the  components  reaches  its  yield  limit.  As  a  strength  design  cri¬ 
terion,  the  PEM  should  not  be  susceptible  to  plastic  deformation 
because  of  its  important  function  [5];  i.e., 


Fintemal-stack  —  Y^PEM  ^rib  5 


(24) 


where  y  is  the  number  of  ribs  on  a  single  side  of  the  BPP.  Combin¬ 
ing  Eq.  (16)  with  (24),  the  maximum  allowable  tightening  torque 
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applied  to  a  single  bolt  when  considering  the  strength  of  the 
internal-stack  is 


T  < 


(^internal-stack  +  ^external-stack  ^V^PEM  ^rib^t 


eq 


^internal -stack 


(25) 


where  cell-core  represents  the  assembled  component  comprising 
BPP-rib,  GDL  and  PEM,  which  were  explained  in  detail  in  our  earlier 
work  [8]. 

3.5.  Strength  design  of  the  clamping  bolts 


3.4.  Thermal  effects 


As  one  of  the  design  requirements  is  that  no  component  gives 
rise  to  plastic  deformation  in  any  operation  situation  (in  other 
words,  only  elastic  deformations  occur),  the  influence  of  the  ther¬ 
mally  induced  deformation  does  not  depend  on  the  history  of  the 
temperature  variation,  but  merely  depends  on  the  magnitude  of 
the  temperature  variation  from  the  assembly  to  the  stable  operat¬ 
ing  condition.  When  taking  washers  into  consideration  (see  Fig.  6), 
our  original  equations  [8]  for  calculating  thermal  effects  require 
amending  as  follows: 


F internal-stack  ~  / 


Finternal-stackFclamping 


Nnternal-stack 


+  k( 


external-stack 


,  /,  r  y thermal  _( ythermal  ,  y thermal 

“I- 1  ML  ^whole-stack  ^internal-stack  '  ^washer  'F 


(26) 


n  _  ^extemaZ-stack^c/amping 

F external-stack  —  T  ~T 

K internal-stack  '  Kextemal-stack 

i  i.  \y thermal  _r  ythermal  ,  y thermal \  1  (971 

'  B  |_  ■ whole-stack  ^ external-stack  '  z washer  ']  ’  ^  ' 

where  positive  values  of  Fintemai_stack  and  Fextemal_stack  represent 
tensile  loads,  zthemml  represents  thermal  deformation  with  the  sub- 
script  denoting  the  corresponding  component,  and  kA  and  kB  are 
two  parameters  merely  used  for  equation  simplification  and  can 
be  written  as 


kA  = 


_ Ckwasber  kjnternal-stack _ 

^washer  +  (kintemal-stack  T  kexternai_stack) 


(28) 


A  clamping  bolt  of  a  PEMFC  stack  mainly  bears  axial  tensile 
stress  that  increases  with  the  tightening  torque.  On  the  basis  of  the 
mechanics  of  materials,  the  axial  tensile  stress  in  each  bolt  should 
meet  the  design  criterion  of  the  material  strength: 


® bolt  ~ 


4F damping 


Cnd\ 


yield 

b0,t  =Wbo»], 


H-safe 


(34) 


where  d\  is  the  minor  diameter  of  the  bolt,  nsafe  is  the  safety  factor, 
and  [crboit]  is  the  allowable  stress. 


3.6.  Design  based  on  the  structural  strength 

In  the  following,  a  design  demonstration  is  given  for  the  PEMFC 
stack  consisting  of  multiple  single  cells  as  described  in  the  work  by 
Bograchev  et  al.  [5]  (see  Fig.  7).  The  design  numbers  of  cells  are  5, 
10,  20,  30, 40  and  50.  The  required  material  properties  and  opera¬ 
tion  parameters  for  this  demonstration  stack  are  listed  in  Table  1. 
All  fasteners  (bolts,  nuts  and  washers)  are  initially  chosen  as  stan¬ 
dard  series  M8  fasteners.  Thermal  deformations  of  the  washers  are 
neglected  owing  to  their  relative  small  thickness. 

First,  the  equivalent  stiffness  of  the  internal-  and  external-stacks 
with  different  numbers  of  cells  is  calculated  and  shown  in  Fig.  8. 
The  figure  shows  that  for  the  studied  PEMFC  stacks,  the  equiva¬ 
lent  stiffness  of  the  external-stack  is  always  greater  than  that  of 
the  internal-stack,  while  both  decrease  nonlinearly  as  the  design 
number  of  cells  increases. 

As  the  operating  ambient  temperature  changes  frequently  dur¬ 
ing  the  use  of  the  PEMFC  stack,  the  temperature  range  of  each 
component  in  the  stack  also  varies  frequently.  As  modern  PEMFC 


1<b  = 


Ckwasber  kextemal-stack 


Ckwasf ier  +  (kintemal-stack  +  kextemal-stack) 


(29) 


Strictly  speaking,  the  magnitudes  of  the  temperature  variations 
of  different  components  differ.  For  example,  the  temperature  vari¬ 
ation  of  the  internal-stack  ranges  from  the  minimum  environment 
temperature  to  the  maximum  operation  temperature,  while  the 
temperature  variation  range  of  the  external-stack  and  bolts  may 
be  relatively  small  owing  to  natural  heat  transfer.  To  obtain  the 
detailed  temperature  distribution  of  the  whole  stack  requires  a 
numerical  simulation  or  a  temperature  monitor,  which  is  quite 
difficult  and  time  consuming.  However,  the  numerical  simulation 
can  be  replaced  by  a  reasonable  assumption  for  the  tempera¬ 
ture  distribution.  One  of  the  simplest  ways  is  first  to  assume  that 
each  component  has  the  same  temperature  as  the  corresponding 
region  (internal-stack,  external-stack  and  bolts)  and  then  deter¬ 
mine  the  possible  temperature  variations  for  each  region.  With  this 
assumption,  the  equations  for  calculating  thermal  deformations  are 
transformed  as  follows  [8]: 


Ckt,0it 


ythermal 
Z bolt 


I  b ythermal 
'  A.  internal-stack 


ythermal 
' whole-stack 


I  b  ythermal 
'  B  external-stack 


+  fe  +  fcB)Z" 


Ckfroit  +  1<A  +  1<B 


(30) 


ythermal  _  r  n/r  ,  i  \y thermal  ,  n/iythermal 

^ external-stack  ~  '  '  '^BPP -sealant  '  ^ sealant  ’ 


(31) 


ythermal 
^ internal-stack 


—  (M  _L  1  \ythermal  ,  M ythermal 

-  [ivi  l  J^BPP-hase  +  iWZ cell-core  ’ 


(32) 


ythermal  _  yr  ythermal  ,  ythermal  \  ,  ythermal 
* cell-core  ~  BPP-rib  +  ZGDL  )  +  ZP£M 


(33) 


Table  1 

Material  properties  and  operation  parameters  of  the  stacl<  hardware. 


Description 

Value 

Source 

Material  property 

PEM  (Nafion®) 

Young’s  modulus  (MPa) 

190 

[21] 

Coefficient  of  linear  expansion  (I<-1 ) 

90 x 10“6 

[22] 

GDL  (carbon  paper) 

Young’s  modulus  (GPa) 

10 

[21] 

Coefficient  of  linear  expansion  (I<-1 ) 

7.9  x  10“6 

[22] 

Bipolar  plate  (graphite) 

Young’s  modulus  (GPa) 

10 

[21] 

Coefficient  of  linear  expansion  (I<-1 ) 

7.9  x  10“6 

[22] 

Sealant  (VMQ) 

Young’s  modulus  (MPa) 

5500 

[4] 

Coefficient  of  linear  expansion  (I<-1 ) 

77  x 10-6 

[22] 

Clamping  bolt  (steel) 

Young’s  modulus  (GPa) 

209 

[23] 

Coefficient  of  linear  expansion  (K-1 ) 

13 x 10“6 

[22] 

Washer  (aluminum  bronze) 

Young’s  modulus  (GPa) 

110 

[23] 

Operation  parameters 

Gas  temperature  (K) 

353 

[24] 

Gas  pressure— anode  (MPa) 

0.15 

[24] 

Gas  pressure— cathode  (MPa) 

0.15 

[24] 

Assembly  temperature  (K) 

293 

Assumed 

Operating  temperature  (K) 

253-353 

Assumed 

Inlet  relative  humidity— anode 

100% 

[24] 

Inlet  relative  humidity— cathode 

100% 

[24] 
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Fig.  7.  Dimensions  of  the  single  cell  structure  of  the  demonstration  stacks.  (The  washers  are  not  drawn  in  this  figure.) 


stacks  usually  operate  at  temperatures  below  353  K  [25],  the  tem¬ 
perature  of  the  MEA  may  range  from  a  low  ambient  temperature 
(for  example  233  K  in  winter)  to  353 1<  (the  working  state),  whereas 
the  bolt  temperature  may  range  from  the  ambient  temperature 
to  a  temperature  higher  than  the  ambient  temperature  but  lower 
than  the  working  temperature  of  the  MEA.  Table  1  lists  some 
possible  temperature  parameters,  while  Table  2  presents  the  five 
possible  cases  for  the  stack  components.  The  first  and  second 
cases  correspond  to  a  low  ambient  temperature  and  high  ambi¬ 
ent  temperature  respectively  (stack  in  a  rest  state);  case  3  assumes 
that  the  stack  is  working  at  a  low  ambient  temperature;  case  4 
assumes  that  the  stack  is  working  at  room  temperature;  while  case 
5  assumes  that  the  stack  is  working  at  a  high  temperature  (a  limit 
case).  All  the  studied  stacks  were  assembled  at  room  temperature 
(293  K). 


Using  the  parameters  given  above,  the  thermal  effects  can 
be  calculated  using  Eqs.  (26)-(33),  where  the  magnitude  of  the 
temperature  variation  for  any  component  equals  the  difference 
between  its  corresponding  working  temperature  (or  the  rest  state 
if  not  working)  and  assembly  temperature.  Fig.  9  shows  that  all  the 
thermal  loads  vary  nonlinearly  with  the  design  number  of  cells. 
Moreover,  both  the  internal-  and  external-stacks  are  subjected  to 
compression  loads  while  the  bolts  are  subjected  to  tensile  loads 
in  cases  1,  3  and  4,  whereas  the  opposite  is  true  in  the  other  two 
cases. 

In  general,  the  PEMFC  stack  should  be  designed  to  meet  the 
requirements  of  the  structural  strength  for  all  possible  conditions; 
i.e.,  the  five  cases  listed  in  Table  2.  Taking  a  specified  case  for 
example,  once  the  design  number  of  cells  has  been  fixed,  the  cor¬ 
responding  lower  bound  of  the  tightening  torque  is  given  by  Eq. 
(22),  while  the  upper  bound  is  the  minimum  of  tightening  torques 
obtained  from  Eqs.  (23)  and  (25).  Considering  all  five  cases  men¬ 
tioned  above,  the  allowable  range  of  the  tightening  torque  for  the 
studied  stack  is  shown  in  Fig.  10.  That  is,  if  the  tightening  torque  is 
designed  within  the  allowable  range,  no  damage  will  occur. 

The  design  is  not  yet  finished  since  the  strength  of  the  clamping 
bolts  has  not  been  checked.  Normally,  only  the  upper  bound  of  the 
tightening  torque  needs  to  be  calculated  using  Eq.  (34)  for  the  bolt 


Table  2 

Assumed  temperature  variations  in  four  cases. 


Case  1 

Case  2 

Case  3 

Case  4 

Case  5 

Assembly  temperature  (K) 

293 

293 

293 

293 

293 

Ambient  temperature  (K) 

233 

313 

233 

293 

353 

Assumed  working  temperature3  (K) 

Internal-stack 

233 

313 

353 

353 

353 

External-stack 

233 

313 

313 

333 

353 

Bolt 

233 

313 

298 

313 

353 

a  Working  temperature  includes  the  temperature  when  the  stack  is  at  rest  in  the 
non-operating  case. 
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Number  of  Cells 


—  □  —  Internal-stack 
— ^ —  External-stack 
— A —  Bolt 


Fig.  9.  Calculated  thermal  loads  for  the  cases  given  in  Table  2. 


strength  design.  For  the  studied  stacks,  it  is  found  that  the  tight¬ 
ening  torque  is  located  in  the  allowable  area  (see  Fig.  10)  and  the 
calculated  safety  factor  for  the  bolt  strength  is  2.8  when  the  bolt 
diameter  is  designed  as  8  mm  (M8). 

4.  Design  based  on  the  stack  performance 

The  optimal  clamping  load  can  be  chosen  in  the  area  of  the 
allowable  range  after  the  strength-based  design.  It  is  an  interest¬ 
ing  and  important  topic  how  to  design  the  optimal  clamping  load 
via  controlling  the  tightening  torque  for  a  special  PEMFC  stack. 
Strictly  speaking,  searching  for  the  optimal  clamping  load  is  a 
multi-objective  optimization  problem,  which  is  a  difficult  topic  in 
the  computational  mechanics  since  the  clamping  load  dramatically 
affects  the  stack  performance  in  a  complex  manner.  One  needs  first 
to  accurately  know  the  relations  between  the  performance  param¬ 
eters  and  the  clamping  load.  Flowever,  on  the  basis  of  the  equivalent 


stiffness  model,  we  introduce  here  a  method  for  obtaining  the  opti¬ 
mal  clamping  load  for  a  given  PEMFC  stack. 

4.1.  Variation  in  the  interfacial  resistance  between  the  GDL  and 
BPP  with  tightening  torque 

The  contact  resistance  of  a  pair  of  solid  surfaces  is  governed 
by  the  multiscale  surface  topography  and  the  surface  physical 
properties.  When  two  solid  surfaces  are  in  contact  under  a  small 
compression  force,  only  the  tips  of  the  surface  asperities  come  into 
contact  and  thus  the  real  contact  area  is  very  small  [26].  In  such 
a  case,  the  interfacial  resistance  is  very  high.  The  real  contact  area 
increases  with  the  compression  pressure,  while  the  Ohmic  resis¬ 
tance  decreases.  For  the  fuel  cells,  the  clamping  load  should  be  large 
enough  to  give  an  acceptably  small  overall  Ohmic  resistance. 

Several  experiment  datasets  for  the  interfacial  contact  resis¬ 
tance  of  a  single  fuel  cell  and  the  total  resistance  including  the 
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Case  1 

Case  2 

Case  3 

Case  4 

Case  5 

Lower  Bound 

- □ - 

—  o  — 

—  A — 

— v — 

— o— 

Upper  Bound 

Fig.  10.  The  allowable  range  of  the  tightening  torque  applied  to  a  single  nut  when 
considering  all  five  cases  listed  in  Table  2.  (Hatched  areas  are  those  for  which  a 
tightening  torque  is  not  allowed.) 


bulk  resistance  are  available  in  the  literature.  Mishra  et  al.  [27] 
measured  the  total  resistance  and  the  interfacial  contact  resistance 
between  different  GDLs  and  graphite  BPPs.  They  proposed  a  fractal 
asperity-based  model  to  predict  the  interfacial  contact  resistance 
as  a  function  of  pressure,  material  properties  and  surface  geom¬ 
etry.  Wang  et  al.  [28]  reported  the  interfacial  contact  resistances 
between  a  Toray  carbon  paper  (Electro-Chem  Inc.)  and  four  types 
of  stainless  steel  BPPs.  Ihonen  et  al.  [29]  conducted  in  situ  mea¬ 
surements  of  the  interfacial  contact  resistance  as  a  function  of 
time,  clamping  pressure,  gas  pressure  and  current  density.  From  the 
experimental  data  given  in  the  work  by  Mishra  et  al.  [27],  Zhou  et 
al.  [1  ]  established  a  mathematical  model  for  the  interfacial  contact 
resistivity  between  a  GDL  and  BPP: 

where  p  is  the  interfacial  contact  resistivity,  ,  B2  and  B3  are  coeffi¬ 
cients  determined  by  experiments,  and  Pcdl  is  the  contact  pressure 
at  the  GDL,  which  can  be  calculated  using 


Pgdl 


^internal-stack  Pgas^GDL 


VAib 


(36) 


Therefore,  the  interfacial  contact  resistivity  can  be  determined  by 
combining  Eqs.  (16),  (35)  and  (36): 


B2yA  rib  deq  (  ^internal-stack  +  ^external-stack  ) 


TCkinternal-stack  PgasAcDL  deq(  k^gr^^siac^  +  k external-stack ) 


(37) 


4.2.  Variation  in  the  GDL’s  in-plane  permeability  with  tightening 
torque 

Permeability  is  another  important  performance  parameter  of 
the  GDL  that  depends  on  the  clamping  load.  In  a  PEMFC,  liquid  water 
is  produced  during  operation  and  exchanges  at  the  interface  of  the 
catalyst  layer  and  microporous  layer.  During  the  water  transporta¬ 
tion,  the  capillary  pressure  gradient  must  overcome  the  negative 
pressure  gradient  of  the  gas  phase.  Therefore,  a  low  permeability 
leads  to  a  high  gas  pressure  gradient  and  high  liquid  water  satura¬ 
tion,  which  in  turn  reduces  the  stack  performance  [30].  Chang  et 
al.  [24]  measured  the  effects  of  the  clamping  load  on  the  perfor¬ 
mance  of  the  PEMFC  and  provided  an  empirical  relation  between 
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Fig.  11.  Variation  of  the  GDL  interfacial  resistivity  and  permeability  with  the  tight¬ 
ening  torque.  (Hatched  areas  are  those  for  which  a  tightening  torque  is  not  allowed; 
the  design  number  of  cells  is  20.) 


the  in-plane  permeability  k  and  the  clamping  load: 

k  =  kq  exp[-8.23  x  103 Pgdl],  (38) 

where  /c0  is  the  initial  in-plane  permeability.  In  the  present  paper, 
Eqs.  ( 1 6),  (36)  and  (38)  are  combined  to  obtain  the  relation  between 
the  GDL  permeability  and  tightening  torque: 


k  —  kq exp 


—8.23  x  10  [TCkinternai  -stack  PgasAcDL  deq  (  ^internal- stack  H^external-stack  )  ] 
yArib  deq  (  k internal-stack H<  external- stack  ) 


} 


(39) 


Taking  a  stack  consisting  of  20  cells  as  an  example,  variations 
in  the  GDL  interfacial  resistivity  and  permeability  with  the  tight¬ 
ening  torque  are  shown  in  Fig.  11.  It  is  shown  that  both  of  the 
parameters  decrease  nonlinearly  with  an  increasing  clamping  load. 
However,  they  are  contradictory  factors  influencing  the  stack  per¬ 
formance.  Thus,  the  optimal  tightening  torque  should  be  chosen 
after  balancing  all  possible  effects.  For  the  present  design  example, 
as  the  interfacial  resistivity  changes  only  slightly  within  the  allow¬ 
able  range  of  the  tightening  torque,  the  optimal  tightening  torque 
can  be  determined  mainly  on  the  basis  of  GDL  permeability;  i.e., 
the  tightening  torque  should  be  applied  as  close  to  the  left  side  of 
the  allowable  region  as  possible. 


5.  Discussions 


5.1.  Considering  the  equivalent  stiffness  of  nuts 

In  general,  the  effects  of  nuts  on  the  equivalent  stiffness  model 
can  be  ignored  for  several  reasons:  (a)  compared  with  the  clamping 
bolts,  especially  in  a  large  PEMFC  stack,  deformations  of  the  nuts 
are  rather  limited,  (b)  unlike  spring  washers,  the  equivalent  stiff¬ 
ness  of  the  nut  is  very  low,  and  (c)  the  clamping  force  mainly  acts 
on  the  first  few  threads  (close  to  the  EP)  [31  ].  However,  when  the 
stiffness  of  the  nut  has  to  be  considered,  we  can  simply  extend  the 
effective  length  of  the  clamping  bolt  by  a  suitable  percentage  of  the 
nut  thickness. 

5.2.  How  thermal  stress  affects  the  optimal  clamping  load 

If  all  thermal  effects  are  neglected,  the  design  process  is  sim¬ 
ple  since  no  working  condition  has  to  be  considered.  However,  the 
thermal  calculation  cannot  be  skipped  in  the  practical  design  owing 
to  its  great  effect  on  the  design  results.  Fig.  12  compares  the  design 
that  considers  the  thermal  effect  and  the  design  that  does  not. 

It  is  found  that  the  allowable  range  of  the  tightening  torque 
when  considering  thermal  effects  is  much  smaller  than  that  when 
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Fig.  12.  Comparison  of  the  design  results  with  and  without  consideration  of  thermal 
effects. 

neglecting  thermal  effects.  In  other  words,  thermal  loads  (stresses) 
are  so  large  that  they  may  damage  the  stack  components  if  the 
clamping  load  is  not  optimized  by  considering  the  thermal  effect. 

5.3.  The  role  of  spring  washer  stiffness 


of  the  washer  may  also  cause  the  stiffness  of  the  whole  stack  sys¬ 
tem  to  decrease  to  an  unacceptable  level  and  thus  may  reduce  the 
vibration  stability  of  the  stack  system.  Therefore,  similar  to  the  case 
for  the  design  of  the  optimal  tightening  torque,  the  optimal  equiv¬ 
alent  stiffness  of  the  washers  has  to  be  determined  by  considering 
all  the  practical  requirements  of  the  stack  system. 

6.  Conclusions 

In  the  present  paper,  a  complete  design  example  of  a  large 
PEMFC  stack  was  introduced  to  explain  step  by  step  how  to  design 
the  optimal  tightening  torque  using  the  equivalent  stiffness  model 
proposed  by  the  present  authors  in  our  previous  work  [8].  The 
design  considers  the  stack  structure  strength,  stack  performance 
and  thermal  effects.  It  was  shown  that  the  equivalent  stiffness 
of  the  internal-stack  and  that  of  the  external-stack  decrease  with 
an  increasing  number  of  cells  in  the  stack  in  a  nonlinear  manner 
and  that  thermal  effects  cannot  be  ignored  in  the  structure  design. 
Moreover,  how  to  obtain  the  optimal  clamping  load  for  a  given 
PEMFC  stack  was  explained  according  to  the  contact  resistance  and 
permeability  of  the  GDL,  and  in  a  practical  design,  the  optimal  tight¬ 
ening  torque  can  be  determined  after  balancing  all  possible  effects. 
Finally,  the  effects  of  the  equivalent  stiffness  of  the  spring  washer 
on  the  structure  stress  in  the  stack  were  discussed.  The  basic  design 
method  described  in  the  present  paper  can  be  used  for  any  PEM 
stack  design  and  may  be  helpful  for  other  fuel  cell  stack  designs. 


Besides  preventing  the  nuts  from  loosening,  the  spring  washers 
(hereafter  simply  referred  to  as  washers)  have  other  useful  func¬ 
tions.  One  function  is  the  adjustment  of  the  thermal  loads  inside 
the  stacks  through  the  design  of  a  suitable  equivalent  stiffness  for 
the  washer. 

Taking  the  stack  consisting  of  20  cells  in  case  4  as  an  example,  the 
effects  of  the  equivalent  stiffness  of  washers  on  the  thermal  loads 
are  shown  in  Fig.  13.  The  figure  shows  that  the  absolute  values  of 
the  thermal  loads  in  both  the  internal-  and  external-stacks  decrease 
nonlinearly  with  an  increase  in  the  equivalent  stiffness  of  the  spring 
washer.  If  the  stiffness  of  the  washer  is  designed  to  be  low  enough 
(less  than  2.5  x  108  N  m-1 ),  the  effects  of  the  thermal  loads  can  be 
effectively  reduced.  However,  if  the  stiffness  exceeds  5  x  1 08  N  m-1 , 
the  washer  has  almost  no  effect  on  the  thermal  loads  acting  on  the 
bolts,  internal-stack  and  external-stack. 

Although  low  stiffness  of  the  washer  is  helpful  in  reducing  the 
thermal  stress  in  the  stack  structure,  an  unreasonably  low  stiffness 
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Fig.  13.  Effects  of  the  equivalent  stiffness  of  the  washer  on  the  thermal  loads  (case 
4;  the  design  number  of  cells  is  20). 
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